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The dynamics after a quantum quench is determined by the weights of the initial state in the 
eigenspectrum of the final Hamiltonian, i.e., by the distribution of overlaps in the energy spectrum. 

We present an analysis of such overlap distributions for quenches of the anisotropy parameter in the 
one-dimensional anisotropic spin-1/2 Heisenberg model (XXZ chain). We provide an overview of the 
form of the overlap distribution for quenches from various initial anisotropies to various final ones, 
using numerical exact diagonalization. We show that if the system is prepared in the antiferromag¬ 
netic Neel state (infinite anisotropy) and released into a non-interacting setup (zero anisotropy, XX 
point) only a small fraction of the final eigenstates gives contributions to the post-quench dynamics, 
and that these eigenstates have identical overlap magnitudes. We derive expressions for the overlaps, 
and present the selection rules that determine the final eigenstates having nonzero overlap. We use 
these results to derive concise expressions for time-dependent quantities (Loschmidt echo, longitu¬ 
dinal and transverse correlators) after the quench. We use perturbative analyses to understand the 
overlap distribution for quenches from infinite to small nonzero anisotropies, and for quenches from 
large to zero anisotropy. 


I. INTRODUCTION AND MOTIVATION 

The past decade has witnessed a tremendous resurgence of interest in the physics of quantum many-body systems 
out of equilibrium [1-3]. This growth has been partly motivated by remarkable developments in experiments, especially 
in experiments on cold atoms, which made the explicit observation of time evolution of many-body quantum systems 
possible [4-14]. In addition to the experimental motivation, fundamental conceptual issues, such as thermalization in 
isolated systems and the nature of adiabaticity, have also played a strong role in driving the development of this field. 
A central paradigm in the study of non-equilibrium physics is the quantum quench. A quantum quench involves a 
protocol starting from the ground state (or a thermal state or another eigenstate) of an ‘initial’ Hamiltonian Hi, and 
then rapidly changing the Hamiltonian so that nontrivial time evolution occurs under a different (‘final’) Hamiltonian 
Hf. Many aspects of the dynamics induced by quantum quenches have been considered in recent years; the literature 
is rapidly growing and already too vast to review here. 

A quantum quench can lead to nontrivial dynamics as the initial state generally has overlap with many eigenstates of 
the final Hamiltonian. If jn) and are the eigenstates and eigenvalues of the final Hamiltonian, i.e., Hf\n) = En\n), 
an initial state [dto) will lead to the following time dependence of the system wavefunction: 

\A>{t)) = ^ Cn\n) = ^ e“'^"‘c„ln) with c„ = (njT'o) . (1) 
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If a large number of final eigenstates have non-negligible overlap jc„j, this sum can lead to highly nontrivial time 
dependence. Clearly, both the overlaps jc„j = [(dto 1*^)1 and the corresponding eigenenergies are important in 
determining the snbsequent dynamics. The object of study of the present paper is the distribution of overlaps, i.e., 
the overlap magnitudes jc„j as a function of energy E„. We will refer to this as the “overlap distribution”. Note that, 
if the squares \cn\^ were considered, this could be regarded as a discrete probability distribution, as J2n ^ and 

the squares jc„j^ represent probabilities. Since jc„j and jc„j^ provide equivalent information, we will focus on jc„j. 
Also note that the distribution of phases of c„ may provide useful information but are not individually well-defined, 
since the phases of individual eigenstates jn) can be chosen arbitrarily. 

This work presents an investigation of the features of the overlap distribution for various quenches in a particnlar 
system. In this introductory section, we provide motivation for the study of the overlap distribution, by reviewing its 
connections to topics that have attracted widespread attention in recent years. 

The overlap distribution is closely related to the statistics of the work performed in a quench, which is characterized 
by the so-called work distribution [15-18]: 

p{w) = 27r^ \cn\'^5{w - {En - Ei)) , 
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where Ei is the initial energy. The work distribution is thus simply a continuous version of the distribution of squared 
overlaps, with the zero of the energy shifted by Ei. The work statistics in non-equilibrium processes in various quantum 
many-body systems is a topic of much recent interest [18-37]. Another quantity of wide current interest, closely related 
to the overlap distribution and work distribution, is the return probability, 

C{t) = |(vI/o|e-*"^*|d/o)P = • (3) 


This is a special case of the Loschmidt echo where the reverse time evolution occurs under the Hamiltonian Eli, for 
which the initial state [dto) is an eigenstate. Following common practice, we will refer to this quantity simply as the 
Loschmidt echo. The behavior of the Loschmidt echo after a quantum quench has been studied in a large body of 
recent literature, see, e.g., [18, 22, 25, 34, 38-76]. In particular, the Loschmidt echo is central to the study of so-called 
dynamical phase transitions [67-76]. By expanding |4'o) in terms of the eigenstates |n) of the final Hamiltonian Hj, 
we find 
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Knowledge of the overlap distribution, i.e., the |c„|’s and corresponding EnS, thus leads directly to the Loschmidt 
echo. For example, the energy width of the distribution determines the coefficient of the initial quadratic decay of 
C{t), as it can be seen by expanding in Eq. (3) to second order in t [46, 48]. 

After a quench, a generic observable O has the following time evolution: 


{0{t)) = (vE-(t)|0|vI/(t)) = (vI/o|e'^^‘Oe-‘^/‘|vI/o) = ^ ^ |c„|2(n|0|n). (5) 
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The long-time behavior of such expectation values, limt_,.oo((!l(<)), is of fundamental interest since this is closely 
related to the question of whether an isolated quantum system thermalizes in some sense. As outlined above, this 
long-time limit is determined solely by the overlap magnitudes and the diagonal matrix elements. This is often called 
the “diagonal ensemble” value of the long-time limit [77]. The eigenstate thermalization hypothesis (ETH), which 
is expected to hold for non-integrable systems, postulates that the diagonal matrix elements (n|C>|n) are smooth 
functions of eigenenergy, for large enough systems [77-81]. The question of thermalization can then reduce to the 
question of how narrow in energy the overlap distribution is. Recent work (such as the “quench action” approach [82- 
88] and an approach based on linked-cluster expansions [89, 90]) has formulated calculations of long-time expectation 
values directly in the thermodynamic limit, thus avoiding the explicit calculation of overlaps in finite-size systems. 
The diagonal ensemble formulation of Eq. (5) for finite sizes, involving the overlap distribution |c„|, will play an 
essential role in the discussion of such calculations for long-time stationary state values of observables, in particular 
for comparing the ordering of the limiting procedures (limits of large sizes and large times). We also note that the 
overlaps c„ appear not only in the long-time limit, but also in the arbitrary-time expression for {0{t)) (cf. Eq. (5)), 
although in this case the phases of c„ are also important. 

From the discussion above, it is clear that the overlap distribution plays a central role in many aspects of quench 
dynamics. Numerical data on overlap distributions have appeared in many different non-equilibrium studies for specific 
quenches, e.g., in Refs. [41-47] where the overlap distributions are called “local density of states”, in Refs. [77, 91-96] 
in the context of thermalization, in Ref. [20] in connection with intermediate-time behaviors of the Loschmidt echo, 
and in Ref. [97] for a spin spiral initial state. There have been, however, relatively few systematic studies of the overlap 
distribution so far. The present work is a step towards addressing this gap in the literature. 

We provide a detailed study of the overlap distribution for a well-known model system, namely the one-dimensional 
anisotropic spin-1/2 Heisenberg model (XXZ chain) at zero magnetization. We consider quenches of the anisotropy 
parameter A, see Eq. (6), and provide an overview of the overlap distributions obtained in quenches from various 
initial values to various final Ay values, using full numerical diagonalization. For the special case of A^ = oo and 
Ay = 0, i.e., quenching from the Neel state to the so-called XX point, we present a full analytical study of the overlap 
distribution. The XX point is mappable to a chain of free fermions, which allows this detailed analysis. We find that 
all nonzero overlaps are exactly equal and derive expressions for them, using the Slater determinant structure of 
the XX eigenfunctions, as well as using the Bethe ansatz. We derive the selection rules determining the eigenstates 
that have nonzero overlaps for this quench, both in the language of fermionic momenta and using the language of 
Bethe roots. We show how these results on the overlap distribution can be used to derive explicit expressions for 
the Loschmidt echo C{t) and for equal-time correlators after the quench. For quenches ending at small but nonzero 
anisotropy. Ay ^ 1, a splitting structure appears in the overlap distribution; we analyze the main features of this 
structure using perturbation theory. Similarly, for quenches to the XX point starting from large but finite initial value 
Ai ^ 1, a perturbative expansion in l/A^ explains the structure of the overlap distribution. 
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After introducing the model and some notation in Sec. II, we provide our numerical overview of overlap distributions 
in quenches from various to various A/ in Sec. III. Section IV is devoted to a detailed analysis of the A^ = oo —)■ 
Af = 0 case. We also use the results for the overlap distribution to derive analytical expressions for time-dependent 
quantities, in particular, for the Loschmidt echo and for longitudinal and transverse two-point correlators. Section V 
outlines a perturbative treatment of the Ai = oo—?>A/<Cl case (perturbation in Af), while Sec. VI describes a 
perturbative treatment of the Ai3>I—?>A/ = 0 case (perturbation in l/A^). 


II. THE MODEL 

This work focuses on the one-dimensional anisotropic spin-1/2 Heisenberg model (XXZ chain) with nearest-neighbor 
interactions. The Hamiltonian is given by 

L L L 

H = Y. = 2 E (E^7+i + ^ E ’ (6) 

i=i i=i 

where S'" = (t“/ 2, a = a;, y, z, with cr" being the usual Pauli matrices, and Sj*^ = SJ ± iSj. We will refer to the A = 0 

part of the Hamiltonian as the XX Hamiltonian: Hxx = ^j+i ■ Throughout this manuscript 

we choose the number L of lattice sites to be even. We also distinguish between two different types of boundary 
conditions. For an open chain we set = 0, for a periodic chain = Sf, a = x,y, z. 

The XXZ chain is a natural generalization of the isotropic Heisenberg chain [98] . It is one of the most widely studied 
models in quantum magnetism in particular, and condensed matter physics in general. It is absolutely central to the 
study of integrability, having a relatively simple integrable structure, yet possessing rich physics and displaying a 
phase transition. In the current era of non-equilibrium physics, it has been widely used as a model system to explore 
and exemplify non-equilibrium phenomena, including studies of the Loschmidt echo (e.g. [22, 46, 64-68]), studies of 
quenches of the anisotropy parameter A (e.g. [99, 100]), and of quenches starting from a Neel state (e.g. [46, 68, 84- 
86, 99, 101]). Dynamics starting from various non-uniform states (such as domain wall initial states), leading to 
propagation or transmission along the XXZ chain, is also a rapidly growing topic of investigation (see, e.g., [22, 102- 
117]). The XXZ chain is also being increasingly used to study real-time thermal transport and real-time dynamics at 
finite temperatures [113-118]. It is a primary testbed for developments and applications of Bethe ansatz techniques 
for addressing non-equilibrium issues (e.g. [22, 82-86, 102, 106, 119-124]). We therefore find the XXZ chain to be an 
appropriate model for a systematic study of the work distribution in quantum quenches. 

In this work, we restrict to positive anisotropies A > 0. The isotropic critical point A = 1 separates two regimes 
with different behaviors; in terms of low-energy physics, the A < 1 regime is gapless, while the A > 1 regime is 
gapped. The Hamiltonian (6) commutes with the z-component total spin operator. For even 

L the zero magnetization sector contains the ground state for all A > 1. Our initial states are ground states of the 
XXZ chain and, hence, lie in the zero magnetization sector. Thus, we can restrict to the subspace spanned by energy 
eigenstates jn) with S^\n) = 0. At A = 0 the chain can be mapped onto a gas of free fermions via the Jordan-Wigner 
mapping, and subsequently solved using free-fermion techniques. For generic values of the anisotropy parameter A the 
XXZ model is integrable and can be solved by means of Bethe Ansatz techniques [125-129]. In the fermionic language 
the anisotropy A acts as the strength of nearest-neighbor interactions. 

A particular focus of this work will be on quenches starting from the purely antiferromagnetic regime (A^ —)■ oo). 
In the limit A^ —>■ oo the initial Hamiltonian effectively becomes the Ising Hamiltonian, Hi ~ Flising = ^f=i 
The ground state of iJismg is exactly degenerate, with both symmetric and antisymmetric combinations of the Neel 
states having the same energy. For A^ large but finite this degeneracy is lifted. We denote the two Neel states as 

\Ni) = and JiVa) = • (7) 

For Ai = oo, any linear combination of \Ni) and jA 2 ) is a valid ground state and, in principle, any such linear 
combination could be chosen as the initial state with which to calculate overlaps. Physically, we find it most reason¬ 
able to choose the symmetric combination for L/2 even and the antisymmetric combination for L/2 odd, which is 
adiabatically connected to the (unique) ground state for Ai large but finite, 

I'l’a) = ^i\Ni) + a\N 2 )) with a = (-1)^/^ . (8) 

For the periodic chain we set a = e'^® with pq = 0, tt. Since T\^„) = e‘^° Id/o.), where T is the translation operator by 
one lattice site, po is nothing but the momentum of the initial state. For the open chain, a is related to the symmetry 
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of I'I'ct) under reflection, ^|4'o-) = where the reflection operator ^ maps lattice site j to L + 1 — j for all 

Recently, some progress has been made in calculating overlaps using the Bethe ansatz [122-124, 130]. There are 
now exact formulas for the overlap between a Neel state and any eigenstate in terms of the Bethe roots describing 
the eigenstate [122-124, 131, 132]. By taking the thermodynamic limit of the results of Ref. [122], it was possible to 
investigate the long-time limit in quenches to the antiferromagnetic gapped phase A/ > 1 [84-86] and also to the 
isotropic point Aj = 1 [86], using the quench action approach [82]. However, extracting all the overlaps for large but 
finite chains by means of the exact formulas of Refs. [122, 124] is technically quite challenging. For this reason, in this 
work we will use these forumlas only in the limit A/ —> 0, where the structure of Bethe roots is well understood. 


III. FEATURES OF THE OVERLAP DISTRIBUTION: NUMERICAL OVERVIEW 

In this section we present an overview of the features of overlap distributions for quenches from various A^ to 
various Ay, obtained using numerical exact diagonalization. We present overlaps with only those eigenstates of the 
post-quench Hamiltonian Hf that lie in the zero magnetization sector. Since the initial states are always in this sector, 
eigenstates of Hf in other magnetization sectors have zero overlap due to symmetry. 

In cases where eigenstates of Hf are (nearly or exactly) degenerate, numerical diagonalization might give overlaps 
with arbitrary linear combinations of the (nearly) degenerate eigenstates. In all cases where we discuss visible features 
of the overlap distribution, we have taken care to resolve such degeneracies. There are more degeneracies for periodic 
chains, so we show mainly open-chain data. We have checked in each case that the qualitative features are very similar 
for open and periodic chains. 


A. Quench from Ai = oo to A/ = 0 

As explained in Sec. H, for quenches starting from A^ = oo we choose as the initial state j^ko) the combination 
I4>cr) of Neel states as defined in Eq. (8). Overlap distributions for A^ = oo to Ay = 0 quenches are shown in Fig. 1, 
for both open and periodic chains. The remarkable feature of these distributions is that the overlaps are either zero 
or have exactly equal nonzero value. There is no reason to expect that the ground state at A = oo should have any 
particular preference between low-energy and high-energy eigenstates of the A = 0 Hamiltonian. Nevertheless, such a 
physical argument does not sufhce to predict that the nonzero overlaps should all be exactly equal. This is a particular 
feature of the Neel-to-XX quench in the XXZ chain. 

Among the (^^ 2 ) eigenstates in the zero magnetization sector, only have nonzero c„’s. The number of 

nonzero overlaps fixes the nonzero values to be jc„j = 2“^/^+^/^, by virtue of normalization and the fact that they 
are all equal. The counting can sometimes be tricky due to degeneracies. For example, there appears to be only five 
instead of eight nonzero overlaps in the L = 8 periodic case in Fig. 1(c). This is because the symbol at the center 
represents four exactly degenerate eigenstates, which cannot be distinguished. There are also few eigenstates with 
the same eigenenergy but zero overlap. Of course, any linear combination of the degenerate eigenstates is also a valid 
eigenstate. So, there is some ambiguity in the overlaps plotted. By choosing the eigenstates to be the appropriate 
‘Slater determinants’ in fermionic language (see Sec. IV), the overlaps turn out to be either zero or jc„l = 2“^/'*+^/^. 

The overlap distribution for quenches from A^ = 00 to Ay = 0 will be analyzed in detail in Sec. IV, in particular 
we will derive the selection rules, which determine whether or not an eigenstate has nonzero overlap. We will also 
exploit this analysis to obtain closed expressions for time-dependent quantities. 


B. Quenches from Ai = 00 to various nonzero Ay 

We now consider quenches starting from the ['ho-) state (A^ = 00 ) to various nonzero values of Ay. In Fig. 2(a,c,d), 
we show overlap distributions for quenches to Ay values in the gapless region, Ay < 1. When Ay is very small, the 
distribution is only slightly distorted from the Ay = 0 case, as expected from perturbative arguments. However, the 
pattern of distortion is rather remarkable. The 2^/^“^ nonzero overlaps now split into two well-defined branches; one 
branch is shifted upwards and one branch is shifted downwards, see Fig. 2(a). The upward- (resp. downward-) shifted 
branch is biased toward lower (resp. higher) energies. In particular, among the eigenstates satisfying the selection rule, 
the lowest few all have positive shift, the highest few always have negative shift, and the eigenstates with intermediate 
energies can have either positive or negative shift. In addition, many of the overlaps that were zero for the Ay = 0 
case, now have nonzero (albeit tiny) jc„j values. 
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FIG. 1: Overlap distributions for the quench Ai = oo —>■ A/ = 0. Data are shown for both open and periodic chains, for two 
different system sizes. There are nonzero overlaps, all exactly equal to which is 2“®^^ ~ 0.354 for L = 8 

and 2“^ = 0.125 for L = 14. In the periodic cases some of the nonzero overlaps are not visually distinguishable because of 
degeneracies, e.g., in panel (c) there are four degenerate nonzero overlaps at En = 0. 
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We will refer to the difference |c„| — 2 r./4+i/2 absolute value of the Neel overlap with an eigenstate of the 

Ay ^ 1 Hamiltonian and the Neel overlap at Ay = 0, i.e., as the “shift” of the overlap. The splitting 

structure is well-defined as long as the shift is much smaller than otherwise the lower branch is difficult 

to distinguish from the states with near-zero overlap. It is useful to normalize the shift with respect to the Ay = 0 
value, i.e., for each eigenstate connected to one having nonzero overlap at Ay = 0, we define 


normalized shift = 


shift 

2-L/4-H1/2 


|c„|-2-V4+i/2 
- 2-L/A+i/2 


(9) 


In Fig. 2(b), the magnitudes of the normalized shifts are averaged over the eigenstates in each branch. We see this to 
be linear with Ay at hxed L, and linear with L at fixed Ay, i.e., proportional to LAf. This indicates that the shift can 
be calculated perturbatively in Ay for any finite chain. However as the chain size gets larger the perturbative region 
of Ay shrinks. In Sec. V we will analyze the splitting and shift using first-order perturbation theory. The branches 
have internal structures, e.g., the overlap values roughly increase with energy within each branch. 

For a hxed L, as Ay increases further, the splitting structures gets diffuse, and the separation between the 2^/^ 
larger values and the rest smaller values also gets obscured, see Fig. 2(c). Eventually, there is no visible remnant of 
the characteristic structure at Ay = 0, as seen for Ay = 0.5 in Fig. 2(d). For Ay = 0.5 the overlap distribution has 
no striking pattern, and no obvious remnant of the hat distribution at Ay = 0. A slight bias toward smaller energies 
is visible. This may be expected because, for any hnite Ay > 0, the Ai = oo ground state Idtc) should have more 
weight in the lower part of the spectrum than in the upper part. 

This trend gets more prominent at larger Ay. In Fig. 3 we show some examples for Ay > I. The bias toward 
lower-energy states gets progressively stronger. This indicates that for quenches from A^ = oo to Ay > I, the 
dynamical evolution will be determined mostly by low-energy states. For Ay > 3 the ground state dominates the 
overlap distribution very strongly. 


C. Quenches from A; S> 1 to A/ = 0 

We now consider quenches to Ay = 0, starting from A^ ^ I large but hnite. The initial state jlko) is now no longer 
exactly the combination Itko.) of Neel states. From perturbative intuition, we expect that as Ai decreases toward 
smaller values the overlap distribution will deviate gradually from the characteristic shape of Fig. I. This is seen to be 
true in Fig. 4, where we have plotted overlap distributions for an open chain with L — 12, for three different Ai ^ I. 
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FIG. 2: (a,c,d) Overlap distributions for Ai = oo (initial state |4'o) = l^'tr)) in an open chain with L = 14 sites. The splitting 
that can be observed for A/ <C 1 disappears as A/ grows, (b) Magnitude of the average normalized shifts for the two branches 
as function of A/ (double logarithmic scale) and as function of system size L (linear scale). The normalized shift for small A/ 
is observed to be linear in A/ and L. 


The nonzero values exhibit a roughly linear dependence on En as is decreased from co. One expects the 

ground state at any finite to be closer to the low-energy eigenstates of the Ay = 0 than to the higher-energy 
eigenstates. Hence, it is not surprising that the slope of the line is negative. The set of overlaps fall on a very nearly 
linear curve for large A^. For moderate A^, see Fig. 4(c), the overlaps follow a more curved line and show more scatter. 
In addition the near-zero overlap values start to be noticeable. 

The shift of the ground state overlap from can be seen numerically to scale as oc l/A^ for a given L and 

large A^, see Fig. 4(d). The overlap of the highest-energy eigenstate obeying the selection rules scales in the same 
way. Also, the normalized shift for the ground state scales linearly with chain length L for a given A^. 

In Sec. VI, treating the XX part of the Hamiltonian as a perturbation to the Ising Hamiltonian, we will show that 
the shift of the overlaps is exactly a linear function of energy at first order in l/A^. In fact, the normalized shift 
for the eigenstates satisfying the selection rule will be shown to be at first order (exact for periodic chains, 

approximate up to order 1/L for open chains), consistent with the numerical data. 


D. Quenches from Ai = 0 to various A/ 

We now consider quenches starting from A^ = 0. Thus, the initial state is the ground state of a free fermionic 
Hamiltonian. Examples for three Ay values are shown in Fig. 5 for an L = 14 open chain. For any finite value of Ay, 
we expect a bias toward low-energy states. This is seen in all three cases. The bias is naturally stronger for smaller 
Ay . We see that up to Ay = 2 the ground state overlap dominates the overlap distribution. Other than the ground 
state, most of the weight for Ay « 1 goes to a few final eigenstates at relatively low energies. We also see from the 
logarithmic plots in the insets of Fig. 5(a,b) that there is a characteristic decaying shape. The upper envelope of this 
scatter-plot decays roughly exponentially, although a sharp qualitative statement about the decay rate is difficult to 
make, due to finite-size effects. 

The Ay ^ 1 case has some bunching structure reflecting the fact that the energy spectrum near the Ising limit is 
grouped into ‘bands’. Each band is characterized by the number of domain walls in the configurations that dominate the 
eingenstates of that band. Accordingly, the overlap distribution is broken into groups. Interestingly, within each band 
the overlap distribution is dominated by the lower-energy part of that band. Moreover, each band of the distribution 


















7 



FIG. 3: Overlap distributions for Ai = oo (initial state l^'o) = and Af > 1 in an open chain with L = 14 sites. The 

overlap distributions decay with increasing energy. The insets show the same data in a logarithmic scale. 
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FIG. 4: (a,b,c) Overlap distributions for Ai 3> 1 and A/ = 0 in an open chain with L = 12 sites, (d) The normalized shift of 
the lowest-energy eigenstate, as a function Ai for fixed chain length L (double logarithmic scale), and as a function of chain 
length L for fixed Ai (linear scale). 




FIG. 5: Overlap distributions for Ai = 0 and various Af in an open chain with L = 14 sites. In the last panel we see how for 
large values of A/ the energy spectrum is organized in ‘bands’. 




































has roughly similar shape as the characteristic shape of the full distribution for the Af « 1 case (see Figs. 5(a,b)). 
The inset to Fig. 5(c) exemplifies this by zooming into one of the bands. 


IV. ANALYTICAL RESULTS FOR THE QUENCH FROM = oo TO A/ = 0 

In this section, we will present a series of analytical results on the quench from = oo to Ay- = 0, i.e., the quench 
from the combination of Neel states to the XX chain. This is the case addressed numerically in Sec. Ill A and 
Fig. 1. In Sec. IV A, we set notation by providing a brief reminder of how the XX Hamiltonian is diagonlized by 
mapping to Jordan-Wigner fermions. 

In Sec. IV B, we use the fermionic language to derive selection rules for an XX eigenstate to have nonzero overlap 
with |d^cr). The rules are expressed in terms of the momenta (periodic chain) or pseudo-momenta (open chain) of the 
fermions. Physically, one rule arises from the structure of the Neel states |A^i, 2 ), and a second rule arises from the 
well-defined symmetry of the combination l^'cr). We also express the overlap magnitudes as a sum of determinants 
(one for each |A^i, 2 ))- This is evaluated as a function of the chain length L, both by exploiting determinant properties, 
and from normalization, using the fact that all the nonzero overalps are equal for this quench. 

In Sec. IV C, we use the overlap formulas of Refs. [122, 124] (valid for all Af) and take the limit A/ —>■ 0. In this 
way we formulate selection rules in the language of Bethe roots. We show how these rules are equivalent to the rules 
in terms of fermionic (pseudo-)momenta. 

Finally, we provide two applications of our detailed understanding of the overlap distribution in this particular 
case. In Sec. IV D, we derive closed expressions for the Loschmidt echo after the quench. In Sec. IVE, we present 
closed expressions for the full time dependence of equal-time two-point correlators. We consider both longidudinal 
((S'|S'|_i_jj)) and transverse {(S)^ ‘^J+n)) correlators. The expressions for the transverse correlators appear here for the 
first time, to the best of our knowledge. 


A. XX eigenstates in free-fermion langnage 


To diagonalize i/xx, i-C., the A = 0 part of the Hamiltonian (6), one uses the Jordan-Wigner (JW) transformation 

i-i 

= i>j + = Cj exp ' ' 


iTrJ2 




“ 2 


I 


( 10 ) 


Here c] and Cj are usual fermionic operators, satisfying anticommutation relations {Cj , cj} = Sij, {c|, c]} = {Cj, c^-} = 0. 
For periodic boundary conditions, the JW transformation turns F/xx to the free fermionic Hamiltonian = 

+ with the boundary condition = (—l)^c^, where N = (N) is the 

number of fermions in the state. Using a Fourier transformed set of operators. 


- — V 


( 11 ) 


jjPBC 

can be written in diagonal form 


H 


PBC 


= , Cfc = cos{k). 


( 12 ) 


The boundary conditions enforce quantization of the momenta k: 


k € Ka = 


27r(m — a) 


(13) 


m—1 


where a = 1/2 for N even, and a = 0 for N odd. 


For the open chain, the JW transformation leads to This Hamiltonian 


diagonalized be introducing the canonical operators 




IS 


cl = 


Z/ + 1 




J 

■ 


(14) 


i=i 
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leading to 


jjOBC 


eqcjcq , 
9 


eg = cos{q ), 


provided that the numbers q are quantized as 


qGQ 


'I L 

nm I 

^ + lim=l 


Since there is no translational invariance in the open chain, we call q “pseudo-momenta”. 


(15) 


(16) 


B. Determinant formulas and selection rules 


In terms of the Jordan-Wigner fermions introduced above, our initial state (8) can be written as 


' L/2 


LI2 


I n 4-1+^ n 




| 0 ), 


(17) 


Vi=i 


i=i 


where |0) is the vacuum state (all spins down). The eigenstates of the XX Hamiltonian in the zero-magnetization 
{N = L/2 fermion) sector can be written for the periodic (resp. open) chain as 


i{fca}4^i) = ■ ■ ■ 4/2 = 44 • ■ • 4/2 


(18) 


with ka € Ka (resp. qa G Q) for all a = 1,..., L/2, and ka ^ kp (resp. qa ^ qp) for a ^ 13. 

The overlap between an eigenstate of the form (18) and the initial state Idler) from Eq. (17) can now be computed 
using Wick’s theorem. For example, with periodic boundary conditions, we obtain 




1 

TIL 

1 

71 


C 1 C 3 . . . Ci_i4i4 ■ ■ ■ 4/2 1 ^) + '^(0|C2C4 ■ • ■ Clc[i4 ■ • ■ ^Il/2 






;;t 


det 

l<j,a<Ll2 


((0|%-l4 1^^)) + i<4<L/2 |0)) 


(19) 

( 20 ) 


Equation (20) follows from Eq. (19) by separate application of Wick’s theorem on the two correlators. Using definition 
(11) this can be simplified to 


For the open chain we similarly obtain 


1 /1\ 


L/4 r 


det 

l<j,a<L/2 




+ (7 det 

l<i.a<L/2 V / 


•\/2 \L -|- 1 


L/4 r 


det ( sin {qa(2j - 1))) -b cr det ( sin (2g„j)) 

l<i,a<lr/2 V ^ V l<i,a<L/2 V ^ V 


l<j,a<L/2 


( 21 ) 


( 22 ) 


The overlaps (21) and (22) can be zero if both determinants vanish, or if their sum exactly cancels. As discussed 
below, it turns out that both cases can occur. 

We first analyze the formula for the periodic chain. We set a = e‘i’“ and simplify the square bracket in Eq. (21), 


det 

l<j,a<L/2 




+ a det 

l<j,a<L/2 




(23) 


Now, if there is a pair of the form {fc, fc -I- tt} C {^q} 47 determinant in Eq. (23) is zero as two columns of the 
matrix are equal, = e^*^^ for all j. The set Ka of all possible momenta, defined in Eq. (13), decomposes into 

L/2 many of such pairs: 


Ka 



27r(m — a) 
L 


27r(m — a) 
L 



(24) 
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with a = 1/2 for L/2 even and a = 0 for L/2 odd. In order to get a non-vanishing determinant, each pair in this 
decomposition has to be occupied by exactly one momentum. This shows that there are only 2^/^ states for which 
the determinant can be nonzero. If we flip one of the momenta by ± 7 r, i.e., choosing the other member of the pair 
(‘pair flip’), the determinant does not change, again due to Therefore, the determinant has always 

the same non-vanishing value for all of these 2^/"^ states. Note that the overlaps (21) might differ by a factor —1 due 
to the different order of momenta in pair-flipped states. 

The other momentum-dependent factor in Eq. (23) is e“'^“=i “ -|- e‘^“, which is zero if the total momentum 

Safa of the eigenstate \{ka}a=i) and the momentum po of the initial state l^'^) differ by tt (modulo 27r). This 
means that only an even number of pair flips are allowed. Hence, the number of states that have nonzero overlap gets 
reduced by a factor of 2 . 

Thus, there are in total eigenstates |{fca}ffa), all with the same overlap (up to a possible global minus sign 

as dicussed above). Since the initial state |4'cr) is normalized, this common value is 


(^.llfcalff,') = 2 -^/ 4 + 1 / 2 ^ 


(25) 


The derivation above did not require any explicit determinant calculation. We also checked that the result 
can be recovered by directly computing the determinant in Eq. (23). The Vandermonde determinant formula 
deti<j_/<„ = ni<j<i<n ~ ^0 allows to express it as a double product that can be further simpli¬ 
fied. Crucial in the evaluation is the formula - g 2 ip 7 r/n^ _ which can be established by noticing that 

UpZliz - ^ and lim.^i ^ = n. 

The analysis of the overlap formula (22) for the open chain is similar. If there is a pair of the form {q, n — q} G 
{<Za}a=i! both determinants are zero as two columns are proportional to each other: sin( 2 ( 7 r — q)j) = —sin{2qj) or 
sin(( 7 r — q){2j — 1)) = sm{q{2j — 1)), respectively. The set Q of all possible pseudo-momenta, see Eq. (16), decomposes 
into L/2 many of such pairs: 



7rm 

L -f 1 


7rm I 
“ L + 1 J ■ 


(26) 


Again, in order to get non-vanishing determinants, each pair in this decomposition has to be occupied by exactly 
one pseudo-momentum. If we flip one of them by ± 7 r, i.e., doing a pair flip, the second determinant does not change, 
while the first one picks up a minus sign. This shows that there are again 2^/^ states for which the determinants 
are nonzero and that they all have the same absolute value. The two determinants are equal for the state where all 
pseudo-moment a are smaller than tt/2 and also for all states that can be obtained from this state by an even number 
of pair flips. In contrast, for states obtained by an odd number of pair flips they differ by a factor —1. If n denotes the 
number of pair flips, i.e., the number of pseudo-momenta larger than 7r/2, the overlap is proportional to (—1)" + cr, 
which is zero for <7 = 1 and n odd or cr = —1 and n even, and it is ±2 otherwise. This supplemental restriction 
also follows from the reflection symmetry of the initial state. Since the reflection operator 3^ commutes with the 
Hamiltonian, the energy eigenstates are also eigenstates of i.e. ^\{qa\^=i) = 'r\{qa}a=i)- Hence, 

(27) 


and the overlap with I'ko-) can only be nonzero if r = cr. (In the fermionic language the action of on a many particle 
state can be deduced from the formula = (—l)'^c^_|_i_j, where N is the total fermion number. Applying this 

to the Neel state we recover the same rule.) Therefore, as in the periodic case, the number of states that have nonzero 
overlap gets reduced by a factor of 2 due to the reflection symmetry. Again, there are eigenstates with overlap 




2-L/4+1/2 


(28) 


Note that this result is in agreement with the ground-state overlap of Ref. [135], obtained by an explicit determinant 
calculation. All other eigenstates (those with a pair { 9,77 — q} or those with the wrong number of momenta larger 
than tt/2) have zero overlap with 

Let us summarize our findings. We have a simple algorithm that allows us to determine which post-quench eigen¬ 
states (18) have nonzero overlap with the initial state j'kcr) and, therefore, contribute to the time evolution and to 
the non-equilibrium dynamics after the quench. For that to happen two conditions need to be fulfilled: 
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(i) The first condition (‘forbidden pairs’) selects the states for which the determinants are individually nonzero, 
i.e. states with momenta or pseudo-momenta satisfying 

\ka — kp \ ^ TT for all a, /3 G {1,..., L/2} (periodic chain), (29) 

Qa + Qp for all a, /3 G {1,..., L/2} (open chain). (30) 

(ii) The second condition (‘number of pair flips’) selects the states for which the two terms in Eqs. (21) and (22) 

do not cancel. This corresponds to selecting eigenstates with the same symmetry (translation or reflection) as 
the initial state. In the periodic case this restriction can be written as ka = Poj which relates the total 

momentum of the eigenstate to the momentum po of the initial state. For the open chain the restriction is 
imposed by the reflection symmetry of the initial state, S^\^a) = o’|d>cr), via cr = (— 1 )", where n is the number 
of pseudo-moment a larger than 7r/2. 

Note that, if we use \Ni) or IN 2 ) as initial state, only the first rule (i) above applies and there are 2^/"^ nonzero 
overlaps. The second rule (ii) arises because I'hg.) respects the translation/reflection symmetry of the Hamiltonian. 


C. Overlaps derived from a Bethe Ansatz approach 


In this subsection, we treat the overlap distribution for the = 00 to = 0 quench using the Bethe Ansatz 
formalism. Recent work has provided expressions for the overlap between Neel states and Bethe eigenstates at arbitrary 
A/ [I22-I24, 131, 132]. We specialize the formulas of Ref. [124] to the XX limit (A/ = 003 ( 7 ) 0; i-®- 7 '^/2)- We 

will show that the condition of having nonzero overlap with an XX Bethe state can be visualized by adjacent corners 
of certain rectangles in the pattern of possible XX Bethe roots. We further relate this condition to the selection rule 
(29) of single particle momenta. The description here is mostly self-contained, but the interested reader may wish to 
consult Refs. [122, 124] for background and details on related matters. 

For simplicity we focus on the case L/2 even. We start our analysis by presenting the solutions of XX Bethe 
equations with twist and constant inhomogeneity (see, e.g., [126]). The twist deformation is related to twisted boundary 
conditions of the spin chain, = e“^S'f, a = x,y,z. It is necessary to introduce a nonzero twist to resolve issues 

with singular terms in the norm formula for Bethe states in the limit 7 —>■ tt/2. We then use ‘twist-deformed’ solutions 
to calculate overlaps by means of the aforementioned formulas. 

The Bethe equations for anisotropy A = 005 ( 7 ) with twist (/ and constant inhomogeneity e read 


sinh(Aj — e + y) 
sinh(Aj — e — y) 



sinh(Aj — Xji + iy) 
sinh(Aj — Xji — iy) ’ 




(31) 


We consider (/,e G K small. For 7 = 7 r /2 and M = L/2 even, the right hand side simplifies to —e“^, which eventually 
leads to 


tanh(Aj - e) = tan(f (nj - i) - ^) , j = l,...,L/2, (32) 

where rij are integers. For small (/> and e, we have — L/4-I-1 < rij < L/4 for j = 1, ..., Mp, and —L/2-I-1 < rij < —L/4 
or L/4-1-1 < rij < L/2 for j = Mp-|-1,..., L/2. Here, Mp and L/2 —Mp are the numbers of Bethe roots with imaginary 
part zero and 7r/2, respectively. We consider the set {Ai, A 2 , A 3 , A 4 } of four possible Bethe roots, corresponding to the 
numbers n, L/2 + 1 — n, —L/2 n, and —n -|- 1 with 1 < n < L/4. For small (/> we can expand them in the small 
parameter 6 = 2 Lcos(^( 2 n-i)) • Setting A = artanh (tan [f (n — |)]) we obtain 

A 3 = —A -l-y-l-<5-|-e, A 2 = A-l-y-|-(5-l-e, 

^ (33) 

A 4 = — A — (5 -I- e, Ai=A —5-l-e. 

These four approximate solutions can be considered as a deformation of the rectangle with corners (A, A -I- i 7 r/ 2 , —A -I- 
i7r/2,-A). 

In the formulas for the overlaps of Ref. [124] — valid also for the so-called (unnormalized) off-shell Bethe states— 
we first send e —?► 0 and afterwards /> —>■ 0. It can be shown that overlaps of the Neel state with normalized Bethe 

r /o 

states identified by a set of Bethe roots {Xj}^/^^ containing an opposite pair of the form (Ai, A 3 ) and/or (A 2 , A 4 ) are 
of order e. Thus, if there are three or four Bethe roots in one rectangle, there is at least one opposite pair and the 
overlap is zero in the limit e —>■ 0. Therefore, to have non-vanishing overlap, there can be at most two XX Bethe roots 
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in each rectangle. Since there are in total L/2 Bethe roots and LjA rectangles, all rectangles are occupied exactly by 
two Bethe roots sitting at two adjacent corners (‘rectangle condition’). In Fig. 6 we show an example for which the 
rectangle condition is fulfilled. 

The overlap formula for Bethe states [122] consists of two parts, a prefactor and a ratio of two determinants. It 
can be shown that the latter is exactly one in the XX limit, if the rectangle condition is fulfilled, zero otherwise. The 
remaining part is the prefactor that simplifies in the XX limit to 






7=7r/2,0 —>-0 


y^tanhOl 

i=i 


+ 7)tanh(/ij - 7) 


L/4 


2 sinh(2/rj 


= 2 


(34) 


i=i 


Here, the Bethe roots , j = 1,..., L/4, belong to a solution of the untwisted homogeneous Bethe Ansatz equations, 
each of them representing one of the L/4 rectangles. Due to the 7r-periodicity of the tangent function we have 


n„,=iC0t2 [^(n, - i)] = 1, 


and finally obtain 



-I 


7—j-0 1 


0 

0 

0 

2-L/4 


for M ^ L/2, 

if there are more (or less) than two roots in one rectangle, 
if there is a rectangle with roots sitting at opposite corners, 
otherwise. 


(35) 


Here, ‘otherwise’ means that all L/4 rectangles are occupied exactly by two roots sitting at adjacent corners, see 
Fig. 6 for an example. 

The rectangle condition allows for four possibilities for each rectangle and hence, in total, for 4^/^ = 2^/^ many 
states, each contributing to the sum of absolute squares of overlaps with a constant 2“^/^. This can be translated into 
the language of momenta of XX Bethe states. The momentum assigned to a Bethe root, which can be interpreted as 
a single particle momentum, is given by 


Pj =l(Aj) 


with 


p[X) = —iln 


sinh(A + i7r/4) 
sinh(A — i7r/4) 


(36) 


The rule that eigenstates with Bethe roots in opposite corners are forbidden translates into the rule that the corre¬ 
sponding momenta cannot differ by ±7r, which is exactly the selection rule (29). 


Im A 




c 

1 1 

1 1 

1 1 

1 1 
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FIG. 6: Example of a Bethe root configuration (full circles) for a chain of length L = 16. The four ‘rectangles’ are indicated by 
the dashed lines. The ‘rectangle condition’ is fulfilled, leading to a nonzero overlap with the Neel state. Each collection of four 
possible roots with equal |ReA| contains exactly two roots sitting on two adjacent corners of the rectangle. 


D. Application: the Loschmidt echo 

In this section we use our findings for the overlap distribution in the quench Ai = oo^Af=0to compute the 
Loschmidt echo (LE) defined in Eq. (3) . We start from the definition of the Loschmidt echo and exploit the selection 
rules of Sec. IV B to obtain an explicit product formula. For the sake of simplicity we only discuss open chains, but 
the periodic case is very similar. We start by considering the case with only one of the Neel states jfVi) or |fV 2 ), before 
deriving the result for the initial state |4>cr) = (l-A^i) + o' \N 2 ))/V^. Recall that the LE is defined as 

C{t) = \{N,.2\e-^^*\N,^2)\\ 


(37) 
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where the Hamiltonian H of the time evolution operator is given in this subsection by the open chain Hamiltonian of 
the free system. Expanding |fVi, 2 ) in a basis of energy eigenstates |n), we obtain 

(fVn2|e-'^‘|iVi,2) = (38) 

n 

where En is the many-body energy and c„ = (A^i, 2 |n.). As in Sec. IV B, we use the sets of single particle momenta 
{qa}a=i label the eigenstates \n). 

From Secs. IVB and IV C we know that the overlap distribution is flat and that if we take |Vi) or IV 2 ) as initial 
state, only selection rule (i) of Sec. IV B must be fulfilled, namely + g/? ^ tt for all indices a and (3. In terms of the 
single particle energies = cos(g), the previous condition for the momenta can be interpreted in the following way: 

e? + -e?' for all g, g' G {ga}^=i • (39) 

Since the total energy of the state is Edg^j^^^^) = ^9- (3^) oan be rewritten as 

L/2 

(Vi,2|e-‘^‘|iVi,2) = 2-^/2^' n , (40) 

a— 1 

where is a sum over all possible states \{qa}a=-i)^ i-O-, a sum over all subsets {qa}^=i C Q, restricted to the 
states/subsets for which Eq. (39) holds. Since half of the possible single particle energies are positive, while the others 
are their negative counterparts, we can rewrite Eq. (40) as 

(V,,2|e--‘|V,2) = n = n ■ (41) 

In this way condition (39) is automatically satisfied, since it is impossible, when expanding the product, to have in a 
single term the same energies with opposite signs. Therefore, for the open chain the LE corresponding to one of the 
Neel states |Ni^ 2 ) is given by the simple formula 


L/2 


m = n 


cos 


t cos 


irm 

lTT 


(42) 


We note that similar formulas have been derived for periodic chains in Refs. [40, 67]. Our result for the open chain can 
be also obtained by expressing the LE as a determinant and then computing it explicitly. This is however unnecessary. 
As we have seen it can be obtained from simple counting arguments based on the selection rules of Sec. IV B. 

We now consider j'Lcr) as initial state. The LE is still given by a formula similar to Eq. (40), 




L/2 


Yl e 


(43) 


where the prime at the sum now means that both selection rules (i) and (ii) of Sec. IV B must be fulfilled. In addition to 
Eq. (39), the set of pseudo-momenta have to satisfy (—1)" = a = (—1)^/^, where n is the number of pseudo-momenta 
with negative energy. Equation (43) can then be recast as 




£,,>0 


n 

eo>0 


— e 


= cos{€qt) + sin(eqt). 


(44) 


eo>0 


£o>0 


A simple way to prove this is to expand the two products in Eq. (44), and check that all states that do not satisfy 
the selection rules (i) and (ii) cancel. Therefore, the LE corresponding to [dtc) is given by 


C{t) = K^'t^je 


-iHt 


= 


L/2 


COS 

( TTIJI \ 

m—l 



L/2 


-b i^/^ 


Sin 


m—l 


t cos 


L + 1 


(45) 


Both results (42) and (45) are plotted in Fig. 7 for several system sizes L. 
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Initial state: Initial state: |T^) 



FIG. 7: Loschimdt echoes in open chains of sizes L — 8, 14, 60 corresponding to the Neel initial state |Ni, 2 ) (left), see Eq. (42) 
and to the initial state (right), see Eq. (45). For the two smaller sizes, a comparison with Loschmidt echoes numerically 
calculated by means of exact diagonalization shows exact agreement with the derived closed formulas. 


E. Application: time-dependent two-point correlation functions 


We now provide a second application of the results on overlap distributions in Sec. IV B for the quench = oo —)■ 
Af = 0. We will compute the full time-dependence of finite distance two-point correlation functions after such a 
quench. For simplicity we will restrict to the periodic chain and list results which are derived in the large-L limit. We 
provide only a skeleton outline of the derivation, pointing out how the selection rules are used. 

The time evolution occurs under the A = 0 Hamiltonian H^x- The equal-time two-point correlators are defined as 


+n;t) = , (46) 

+n;t) = ■ (47) 


We refer to them as the longitudinal and transverse two-point correlators, respectively. We consider as initial state 
Itko) the four different cases |'I>+) = (|iVi) -|- \N 2 ))/V^, Idt-) = (|Vi) — \N 2 ))/V^ and the two Neel states \Ni) and |fV 2 ) 
themselves. Since the state is an eigenstate of the translation operator T, the two-point correlators only depend 
on the distance between the two points, {j, j + n-,t) = G^^'^ (n; t) = ^ 12^=1 G^J~^ (j,j + n;t). Applying the 

Jordan Wigner mapping (10) we obtain 


Giq»;') = 7i:K<’J+n> = 7i:i 

i=i i=i 


J, 


1 


1 


'Wt 2 J 1 2 


= 7 i:< 


1 


i=i 


(n;t) = jY. 


L ^ 
i=i 


n—1 


11(1-27 


t 'I 


m—1 


j+n 


), 


(48) 

(49) 


where from now on the time dependence is implicit, i.e., the expectation values are taken between the time evolved 
states, (•) = (4'(t)| • |4'(i)) with |4'(t)) = e“'^‘ I'l'o)- In the first line we used that the initial state Idto) hes in the zero 
magnetization sector and therefore ^ — 1/2 for nil times t. In the second line we wrote the ‘Jordan-Wigner 

string’ operators as = 1 — 2c] c^-, which can be written as 












P,Q 

P9^Q 


(50) 


It turns out that the first two terms cancel up to order 1/T, so that the last term dominates in the large L limit. Thus, 
the transverse correlator will involve correlators of products like [nm=i [nm=i , where {pm} and {q-m} are 
sets of n distinct single-particle momenta. As an example, the correlator G'^~ {n; t) is a momentum sum over products 
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of phase factors, which can be rewritten as 




n- 

m—l 


.771=1 






= E m\ 


{k},{k’} 


n< 

.771=1 


H' 

_m=l 


IW> 






tm 

2n 


(51) 


The second line has been obtained using the selection rules of Sec. IVB. Since n momenta of the eigenstates \{ka}ai^i) 
that obey condition (29) are fixed, there are ^ = 2^/^“"“^ states with overlap square 2“^/^+^. This 

gives the factor 2“". The factor (1 + (—1)"')/2 comes from the additional selection rule (‘number of pair flips’) for 
the total momentum of the eigenstates |{A:}) and |{A:'}). Hence, the transverse two-point correlator at points j and 
j + n is only nonzero for n even. Here, A is the mapping matrix from the set {pm} to the set {(?m}, given by an 
element of the symmetric group 5". The factor (—1)1"^^! implies that the result is a determinant. Carrying all factors 
and simplifying, we eventually obtain (t > 0) 


^<4n-m'-l(2t)^ 


GJ (n; t) = - det„ ( (2t)) for n even, 


G+;(n;t) = 0 

Here, we have defined the quantity 


for n odd. 


Xn(2t) = ^ E 


2ii cos{pj)—impj 


Similarly, one obtains for the Neel states themselves 


Gn^ UJ + n;t) = - - -det„ . 


(52) 

(53) 


(54) 


(55) 


In the derivation of the formulas for G^^ some 1/L terms were dropped. However, we have found that formulas (52) 
and (53) for G’jl“(n; t) are exact for n = 2,4 and for any hnite even L. 

The longitudinal correlators are simpler to derive, since they have no Jordan-Wigner strings. The result is 

= + (-l)"Jo^(2f)) (56) 

for all four initial states l^/o) = |'I'±), |iVi^ 2 )- 

In the thermodynamic limit, the summation in Eq. (54) becomes the integral that defines the m-th order Bessel 
function. 


(57) 

Thus, the above formulas for G({l~ and G|,* are also valid in the thermodynamic limit, with the substitution —>■ J^. 
Previously, Ref. [99] presented (connected versions of) longitudinal correlators in the thermodynamic limit. To the best 
of our knowledge, the determinant formulas for the more involved transverse correlators have not appeared previously. 


V. QUENCHES FROM Ai = oo TO A/ < 1: PERTURBATIVE TREATMENT OF THE SPLITTING 

Quenches from the initial state |'ho-) (A^ = oo) to small but nonzero values of A/ display a ‘splitting’ of the overlap 
distribution, as discussed in Sec. HIB and shown in Figs. 2(a,b). In this section, we will use first order perturbation 
theory in A^ to characterize this splitting. 


A. First order perturbation expansion 

For the final Hamiltonian H = i7xx + ^fH^, we consider as a perturbative potential V. Unfortunately, the 

unperturbed Hamiltonian is highly degenerate. In particular, among the states which obey the selection rules, 
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there are many which are degenerate. Hence, a full calculation requires the use of degenerate perturbation theory for 
many of the eigenstates. However, our purpose here is to provide a basic understanding of the splitting and the general 
form of the overlap distribution for small nonzero Af. Hence, it will suffice to concentrate on those states which have 
non-degenerate energy eigenvalues, and use the much simpler non-degenerate perturbation theory. We will show that 
this treatment recovers the main features observed numerically. 

For the sake of simplicity we consider a periodic chain with L/2 even, which means that in the zero magnetization 
sector there is an even number N = L/2 of fermions in the system. All the calculations can be trivially extended to 
the odd particle case. In terms of the Jordan-Wigner fermions of Eq. (10), the perturbative potential V = AjH^ can 
be written as 


^ = ~ 0 (”f+i “ 0 = ^ \ ^ > (58) 

where nj = cjc^ . We can drop the constant term and the terms proportional to the particle number N = 

since they represent just shifts of the energy and are irrelevant for a perturbative expansion. Using again momentum 

operators defined in Eq. (11) the only relevant contribution can be written as 

^ E E . (59) 

Carrying out the sum over the spatial indices we find the condition g-l- s = p + r, which is nothing but the momentum 
conservation. After these considerations we can rewrite the perturbative potential for the periodic chain as 

(60) 


Let us now formally write the expression for the corrected overlap. Given |$) the ‘corrected eigenstate’ with 
limA^-fO l^h) = \{ka}a=i}i til® latter being an eigenstate of the XX chain that obeys the selection rules (‘pair condition’ 
(29) and zero total momentum), we obtain 


= V2 



det 

l<j,a<L/2 



Ek — Ek> 



(61) 


where the states |{A:}) = |{A:a}Ei) = \{^a}a=i) eigenstates of the free fermionic Hamiltonian. Here, 

we only choose states |{A:}) for which the energy value E = J2a=i non-degenerate, such that the energy 
difference in the denominator is never zero. Normalizing with (d^o-Hfc}) we obtain 


('f'.Kfc}) Ek-Ek> 


where we used that the determinant under the sum is either zero or always given by the same absolute value 
I deti<yQ,<L/2(e^'^“^)|. The prime at the sum means that it is restricted to the ‘nonzero overlap states’. The per¬ 
turbative potential (60) connects states |{fc}) and |{fc'}) in which just two momenta are exchanged. Due to the 
selection rule (29) we have that 


\ka-k'J='K, Ikp - k'p\ = TT . 

Therefore, the energy difference in the denominator can be written as 


(63) 


Ek — Ek' = cos{ka) + cos(fc^) — cos{k'^) — cos{k'i^) = 2 (cos(A:q) + cos(fc/ 3 )) . 
Furthermore, for the numerator, we obtain after a short calculation 


(64) 


(65) 


The sign of ({fc}|H|{/c'}) is fixed by the order of momenta of the states |{fc}) and |{fc'}). It exactly cancels the sign 
coming from the determinant that we pulled out of the sum in Eq. (61). The relative minus sign in front of the sum 
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in Eq. (62) is due to the factor e'*^^ in Eq. (60) and comes from the fact that the momenta fco,, kp and k'^ differ 
by ±7r. In summary, we find 

|(d/.|d>)|-|(vl/.|{fc})| ^ sin^(^) 

\{'i>a\{k})\ L ^^^cos{k^) + cos{kf,)- 

a<0 

Follwoing the definition (9) in Sec. IIIB, we call the quantity in Eq. (66) the “normalized shift” and the numerator 
of the left-hand side just the “shift”. 

In Fig. 8 we plot the perturbative results (66) for the nonzero overlaps for A/ = 0.01 and L = 48. (We plot 
the predictions for actual overlaps instead of the normalized shifts.) The perturbative treatment only covers non¬ 
degenerate states. Hence, a quantitative comparison with the numerical data of Fig. 2(a) is not appropriate. However, 
non-degenerate perturbation theory clearly reproduces the main features of the data. Figure 8 shows the overall shape 
to be very similar to that seen in the inset to Fig. 2(a), with two branches each having an overall upward slope. The 
perturbative approach explains the shift being proportional to A/ (by providing a non-vanishing contribution at first 
order in Aj^, see above). In addition, we will show in the subsection below that the result (66) implies that the smallest 
and largest normalized shifts grow linearly with L, so that the normalized average branch shifts grow linearly with L, 
as found numerically, see Fig. 2(b). 

The perturbative prediction in Fig. 8 also shows some substructures, e.g., each branch is father split (recursively) 
into sub-branches. It is unclear whether such substructures might survive when the degenerate states are taken into 
account, or when higher-order corrections are included. In addition, substructures in the XXZ spectrum are often 
heavily determined by boundary conditions [110, 112, 120]. So, they can be expected to be different for periodic and 
open boundary conditions. A detailed study of such substructures is beyond the scope of the current work. 



FIG. 8: Overlap magnitudes for a quench from Ai = oo to small A/, obtained from first-order non-degenerate perturbation 
theory. The specific values plotted here are for A/ = 0.01 and L = 48. Only non-degenerate states are included. The dashed line 
depicts the unperturbed (A/ =0) value q'pg horizontal lines (dots and crosses) indicate the extremal values of 

the shift acoording to Eqs. (67) and (68). 


B. L-dependence of smallest and largest shifts 

In order to see how the normalized shift scales with L we compute the sum in Eq. (66) for large system sizes L. 
First of all, we observe that the shift is minimal in absolute value and opposite for the ground state and the highest 
energy state. For the latter the momenta are uniformly distributed between zero and 7r/2 as well as between 37r/2 and 
27r or equivalently, due to the periodicity of the trigonometric functions, between —7r/2 and 7r/2. So, ka = ‘^{a—112) 
with —L/4 -|- 1 < a < L/4 and the sum over momenta turns in the thermodynamic limit into an integral from —7r/2 
to tt/2 with prefactor L/[2'k). Therefore, the double sum of the right-hand side of Eq. (66) can be written as 

W sin» (t^) 

^os{ka) + cos{kp) 

q</3 


^7r/2 pTT/2 


sin" (^) 


I I /7 \ / cdkdq — 

/_„/2 7-V2 cos(fc)-b cos(q) 2 \4 2tt 


(67) 
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Together with the small prefactor 2A f jLva the perturbative potential (60) this eventually leads to a normalized shift 
of —LA/(l/4—l/(27r)) « —0.09085 LA/ for the highest energy state and, similarly, LA/(l/4—l/(27r)) « 0.09085 LA/ 
for the ground state. This means that as we increase the system size L we have to consider smaller and smaller values 
of A/, i.e. A/ <?; 1/L. Otherwise, we are not in the perturbative regime anymore. To summarize, if we increase the 
value of the final anisotropy from A/ = 0toA/ <?; 1/L, a gap opens in the nontrivial part of the operlap distribution, 
which we call “splitting” and which scales as with prefactor A/(l/2 — I/tt), i.e. twice the scaling factor 

of the shift of the ground state. 

Furthermore, the largest value of the shift can be easily computed, too. Our numerical investigation shows that the 
state that belongs to the largest split is basically given by the set of momenta of the highest energy state where the 
two innermost momenta ±7r/L are replaced by the two momenta ±(7r —tt/L). This leads to an additional contribution 
with respect to the dominating part of the highest energy state, i.e. LA/(l/(27r) — 1/4), of 

2A/ sin"(f-^-¥) _ 2A/ s\v? ^ - ^{2m - 1)) 

L ^ cos(7r — tt/L) + cos(fc/ 3 ) L ^ — cos(7r/L) + cos(^(2to — 1)) 

4LA/ 1 4LA/ 

TT^ mfm — 1) TT^ 

( 68 ) 

The first factor 2 in front of the sum comes from the fact that both innermost momenta are replaced, the second factor 
2 from the fact that the remaining momenta kfj are symmetrically distributed around zero (modulo periodicity of 27r). 
Hence, the largest shift also scales like now with the prefactor A/(4/7r^ + l/(27r) — 1/4) « 0.31444 A/. 

In Fig. 8, the large-L predictions for maximum and minimum values of the shifted overlaps, Eqs. (67) and (68), are 
indicated with horizontal dotted/crossed lines. 
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VI. QUENCHES FROM Ai > 1 TO A/ = 0: ANALYTICAL RESULTS 

In this section we use a perturbative approach to explain one of the features of the overlap distribution that can be 
observed in the quench protocol from A^ large but finite to the free fermion point A/ = 0 . As it can be seen in Fig. 4, 
the nontrivial part of the overlap distribution is rather smooth as a function of and decreases roughly linearly 
with En- The slope increases with decreasing A^. This behavior is in sharp contrast to the splitting structures of the 
Ai = oo—>A/<|;1 case. In the following we will use perturbation theory in l/A^ to show that the normalized shift is 
proportional to at first order (for periodic chains), which explains the smooth behavior of the |c„|’s as a function 
of En- We will also compute the quadratic correction by making use of second order perturbation theory. 

For simplicity we consider a periodic chain. In the case of an open chain the main steps of the following calculation 
can be similarly performed, but there are complications because of the edges. The initial state iTg), which is the 
ground state at large but Hnite A^, is adiabatically connected to the state |4'cr), which we have used as the initial 
state for A^ = oo: 

lim I'I'o) = I'I'ct) with a = (—1)^^^. (69) 

Ai—>00 

In order to make the perturbative calculation mathematically stringent, we divide the initial Hamiltonian (and hence¬ 
forth all corresponding energies) by A^. We consider the second term in 

L 1^1 

i7./A, = ^ + A/ ^ 2 

i=i * 


as a small perturbation. 

T /2 T /2 

In order to analyze the overlaps ({fcajaLil^'o), where |{A:a}a=i) are the eigenstates of the final Hamiltonian with 
A/ = 0, perturbatively, we expand the ground state j'l'o) of Hi in the small parameter l/A^, 

Ivko) = 14/1,°)) + A-i 14-W) + A-2 + ... (71) 

The zeroth-order (unperturbed) wavefunction is = I'I'ct)- 

The energy eigenvalue E„ of j'l'o.) is two-fold degenerate at zeroth order and the corresponding other state is simply 
|4>_o). Since the two states |'I'-cr) and |4>o) (with a fixed) are both eigenstates of the translation operator with 
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opposite eigenvalues —u and cr, we deduce ('I'_cr|'I'o) = 0. This means that |dt_cr) does not appear in any expansion 
of l^'o)- We can therefore apply non-degenerate perturbation theory. Using = |d>^) the higher order corrections 

read 


^0 ) = I] 




(nW|gxx|^.) 
E, - Ui") 


(72) 




|m(0)),|n(0))5^|*±„) 


P J^(^) XT’ 


^ ( TT' _ 2 ^ (TP _ ^7'(0)^2 ’ 






In' 


(0))/|vp±„> ) 


(73) 


where En'^ is the eigenvalue of iJz with eigenstate 

We first compute the matrix elements in the numerator of Eq. (72): 




1^0-) - E 2 (E‘^i-l-1 + E+l) ('^|n(0))GSi + O'<^|n( 0 ))GS 2 ) ’ ('^4) 

V j^i V 

where 5'i^2 are the spaces spanned by the states = {S^+ S~ |fVi_ 2 ), j = i.e., the space 


(0) 


spanned by the spin configurations obtained via a single neighboring spin pair flip from a Neel state. The energies E, 
are all the same, independent of the choice of G Si ^2 and given by En'^ = —Ll(^ 1. Therefore, E^, — En'^ = —1 

for all |n(°)) G S' 1^2 and, henceforth, 




1 

71 


5 E 

|n(0))GSi 



|n(0))GS2 



= -i7xx|^a) . 


(75) 


Thus, the first-order term is obtained by application of 77xx on the initial state (a similar trick was used in Ref. [133]). 

T /O 

Using that the states [{A:}) = \{ka}J^i) are eigenstates of the Hamiltonian iXxx with energy Ek we obtain up to 
first order 


{{m.) A, • 


The second order term is a bit more involved. The first term in Eq. (73) can be computed as 


(76) 


E 

|m(0)).|n(0))/|*±„) 


(mW|gxx|n(°)) (nWjgxxl^g) 


= - E E 

|m(0))^|*±„> |n(0)) 


(m(°)| ffxx V(°))gSi + ^V(°))gS2 

E^ - E^°^ 2 V 2 


(77) 


For each fixed state £ Si U S 2 there are of order many terms for which the flipped spins do not connect to 

each other and for which the energy difference is simply Eg. — Em'^ = —2. In contrast, there are only of order L many 
terms for which the energy difference is not given by this value. We drop these terms because they are subleading in 
L. Since the second and the third term in Eq. (73) are also negligible with respect to the leading contribution, we 
obtain 


l'f'[.'^)=-^(^xx)7'I'.) • 


(78) 


The normalized overlap therefore reads, up to second order in l/A* and by neglecting the 1/L contributions. 


({fc}|d/,) aE2'vAJ 


( 79 ) 
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This provides a second order correction for the nontrivial part of the overlap distribution in the quench protocol from 
Ai 1 to A/ = 0. 

Note that the energy is extensive and bounded by \E}.\ < L/tt. So, if the condition L/tt ^ A^ is fulfilled, we 
expect the perturbation series to converge. The data for open chains from numerical exact diagonalization, see Fig. 4, 
is described reasonably well by the large-L perturbative results derived in this section for periodic chains. 


VII. CONCLUSIONS 

Motivated by the observation that the overlap distribution is important for many aspects of quench dynamics, 
in this work we have considered a prominent paradigm model, the XXZ chain, and undertaken a systematic study 
of the properties of the overlap distribution for quenches of the anisotropy parameter A. Using numerical exact 
diagonalization, we have provided a broad overview of the features of the overlap distribution, for quenches from 
various values A^ to various values of A^. 

In particular, we have thoroughly analyzed the quench from A^ = oo to A^^ = 0, i.e., when the initial state is of the 
form (8) (combination of Neel states) and the post-quench Hamiltonian is described by a free-fermion theory through 
Jordan-Wigner transformation. The overlap distribution in this case has the remarkable feature of being ‘flat’: all 
nonzero overlaps have the same value. Using the free-fermion structure of the post-quench eigenstates, it was possible 
to analyze the overlap distribution analytically, to express the overlaps using determinant formulas, and to explicitly 
calculate their values (see Sec. IV). The states that have nonzero overlap can be chosen using certain ‘selection rules’. 
One rule arises from the spatial structure of the Neel states, while a second rule arises from the symmetry of the 
initial state (8), in terms of translation (periodic boundary conditions) or reflection (open boundary conditions). The 
‘selection rules’ for the periodic chain have been compared with results that we derived from recently obtained Bethe 
Ansatz formulas [122-124]. We found that these rules can be equivalently expressed in terms of adjacent corners of 
rectangles in the pattern of possible XX Bethe roots. We also provided two example applications of our detailed 
understanding of the overlap distribution for the Ai = oo—>Ay=0 quench. Using the selection rules, we were able 
to provide concise expressions for the Loschmidt echo (return probability) as well as equal-time two-point correlators 
after the quench. In particular, time-dependent transverse correlators, (S'^5'“|_„), are somewhat nontrivial due to 
Jordan-Wigner strings. We provide determinant formulas for such correlators, derived by making use of the selection 
rules. 

For quenches from A^ = oo to small Ay > 0, the overlap distribution shows a splitting structure which we analyzed 
perturbatively (perturbation in Ay, see Sec. V). For quenches from large A^ ^ 1 to Ay = 0, the overlap distribution 
acquires a simple dependence on energy which we analyzed using perturbation theory in l/A^, see Sec. VI. 

The present work is related in various ways to recent advances in non-equilibrium physics. Many recent publications 
have provided overlap distributions for specific quenches. This motivates systematic studies like the present one. Using 
the Bethe Ansatz, formulas are now available [122-124] for the overlap between Neel states and Bethe eingestates at 
any A. Since the use of these results requires knowledge of the Bethe roots describing all eigenstates, it is technically 
challenging to explore overlap distributions using this approach. Our numerical overview and detailed explicit analysis 
of specific cases thus provides specific contexts for these general technical results. Our applications to Loschmidt 
echoes and correlators can be compared with existing calculations of the Loschmidt echo [40, 67] and of longitudinal 
correlators [99]). In contrast to these approaches, we have used the structure of the overlap distribution to derive time- 
dependent quantities. This has allowed us to obtain some new expressions. Finally, certain overlaps with low energy 
eigenstates in the gapless phase [Ay] < 1 can be studied using field-theoretical methods. In particular, subleading 
terms in such overlaps [134-136] allow to access universal data about the underlying conformal field theory. Our focus 
here, in contrast, has been on the full spectrum, not just low-energy physics. 

This work opens up many new questions. First of all, the present study is by no means exhaustive. There are 
features of the overlap distribution in the XXZ chain which remain to be analyzed in detail. For example, the physical 
nature of the eigenstates which have large overlap for A^ = 0 to Ay of order one quenches (or vice versa), are 
unclear and might deserve further study. Second, in addition to the applications we have provided for the cleanest 
case (Ai = oo —>■ Ay = 0). It might be possible to exploit the observations on overlap distributions in other quench 
protocols (other A^ and Ay) to extract properties of the time dependence of Loschmidt echoes and correlation 
functions. Finally, it would certainly be interesting to analyze overlap distributions in a series of other models. One 
could ask, for example, whether integrability plays a general role in determining the features of the overlap distribution. 
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